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O '. Abstract 



We consider mathematical aspects of the axion electrodynamics in application to the problem 
of evolution of geomagnetic and terrestrial electric fields, which are coupled by relic axions born 
in the early Universe and (hypothetically) forming now the cold dark matter. We find axionic 
analogs of the Debye potentials, well-known in the standard Faraday - Maxwell electrodynamics, 
and discuss exact solutions to the equations of the axion electrodynamics describing the state of 
axionically coupled electric and magnetic fields in a spherical resonator Earth- Ionosphere. We 
focus on the properties of the specific electric and magnetic oscillations, which appeared as a 
result of the axion-photon coupling in the dark matter environment. We indicate such electric 
and magnetic field configurations as longitudinal electro-magnetic clusters. 
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1 Introduction 



1.1 On the physical aspects of the problem of the axion - photon 
coupling 

The term dark matter axions, which came into general use during the last decade, reflects 
a synthesis of two trends based on ideas of modern astrophysics and cosmology on the one 
hand, and of high-energy particle physics on the other hand. Dark matter, which neither 
emits nor scatters electromagnetic radiation, is considered to be one of the main cosmic 
substrates accumulating about 23% of the Universe energy Mass density distribution of 
the dark matter is now well-studied due to observations and theoretical modeling (see, 
e.g., PQ-I1] for details, review and references). The origin of the dark matter is not yet 
established, but there are several hypotheses about its nature. The most attractive is the 
hypothesis that the dark matter is formed by massive pseudo-Goldstone bosons (axions) 
belonging to the class of WIMPs (Weakly Interacting Massive Particles). The story of 
axions, which were predicted by Peccei and Quinn [5] and introduced into the high-energy 
physics as new light bosons by Weinberg jB] and Wilczek [7], is well-described: one can 
find reviews and references, e.g., in [8]-|15j. 

We are interested in studying of the axion-photon coupling; more precisely, we consider 
the response of the photon system to the action of the axions and possible experimental 
consequences of such action. First of all, it is worth reminding what type of axions we 
are going to consider in this paper. It is well-known that axions can be generated if 
there exist strong magnetic field, strong electric field, and they are not orthogonal to one 
another. Such axion sources can exist in the Solar plasma, in the magnetosphere of neutron 
stars, and in laboratories on the Earth. Other mechanisms of the axion production are 
predicted to occur in the early Universe; the axions born in the early Universe form now 
the unique cosmological substrate indicated as relic (or primordial, or background) axions; 
these are the cosmological dark matter axions, and their density distribution depends on 
time. Since the axions are massive particles, their gravitational attraction forms non- 
homogeneous configurations; for instance, the density distribution in the galactic halos 
is stationary but depends on the distance from the galactic center. In other words we 
live (hypothetically) in an axionic environment and have a right to pose a question: how 
serious the influence of the axions on the terrestrial electrodynamic systems is. 

Wave propagation in the axion electrodynamics is studied in detail for two cases: first, 
when one deals with running plane waves, second, when one uses the approximation of 
geometrical optics, i.e., the problem of dispersion relations can be effectively reduced to the 
analysis of the Fresnel equation (see, e.g., [THl El HB] ) • It was shown that the axion-photon 



interactions induce the effect of polarization rotation, if the electromagnetic waves travel 
through the axion system (see, e.g., [19]- [22] for details). The inertial-gradient extension 
of the Einstein-Maxwell-axion model |23j predicts that the optical activity induced by 
relic cosmological axions should be a nonstationary process and should keep information 
about the rate of the Universe expansion. Nonminimal extension of the Einstein-Maxwell- 
axion model [23] allows us to describe the birefringent wave propagation. Minimal two- 
parameter gradient-type extension of the Einstein-Maxwell-axion model [25J also predicts 
the existence of birefringence, as well as, could explain abnormal growth of the number 
of axions in the early Universe. 

In this paper we focus on the electrodynamic effects, which can occur in the terrestrial 
electric and magnetic fields under the influence of the relic dark matter axions. Relic 
axion distribution is considered to be homogeneous but time dependent. We work here 
within the Weinberg - Wilczek - Ni version of the axion electrodynamics, and focus on 
the axionically coupled electromagnetic oscillations in a spherical resonator. The obtained 
results are new, they can have numerous applications and are planned to be generalized 
in the framework of nonminimal and gradient-type extensions of the Einstein - Maxwell 
- axion model. 

1.2 On the mathematical aspects of the problem of the axion - 
photon coupling 

The axion-photon interaction can be described both in terms of high-energy physics and 
in terms of macroscopic axion electrodynamics. We analyze the problem in the framework 
of axion electrodynamics in the Weinberg - Wilczek - Ni (WWN) version [61 [7J [261 127] . 
The WWN model is the particular case of the Carroll - Field - Jackiw (CFJ) model [28] . 
based on the Chern - Simons ansatz [29]. These models are often discussed in context 
of the Equivalence Principle and Lorentz violation (see, e.g., [SU] - [23] for review and 
references). Two models: CFJ and WWN, are known to be effectively coinciding, when 
the pseudovector pi appearing in the first model is proportional to the four-gradient of 
the pseudoscalar field 0, on which the WWN model is based (e.g., Pi = V«0). Let us 
mention that in the cosmological context, when pseudoscalar field is the function of time 
only, <f>(t), and the pseudovector pi is time-like and has only one component, pi = 5®il>(t), 
both electrodynamic extensions: WWN and CFJ, - give the same results. It gave us a 
supplementary motivation to use the WWN version of the axion electrodynamics in this 
paper. 

Generally, the axion electrodynamics in the Weinberg - Wilczek - Ni version is more 



sophisticated from the mathematical point of view than the Faraday - Maxwell theory. 
The main extension is that the axion electrodynamics includes additional equation for the 
pseudoscalar (axion) field 0, which, in its turn, enters the equations for the electromag- 
netic field strength F^, (the Maxwell tensor), as well as, the constitutive equations. The 
total self-consistent system of master equations is nonlinear. The set of boundary and 
initial conditions is extended accordingly. We are interested to solve the master equa- 
tions of the axion electrodynamics, which guide the evolution of the magnetic and electric 
field in the spherical resonator " Earth- Ionosphere" . In the framework of the Faraday - 
Maxwell electrodynamics this problem is solved and discussed in detail, in particular, in 
terms of Hertz and Debye potentials [35] . These results are classical and they were used 
in numerous applications, in particular, in radio- and tele-communication, navigation, 
magnetic storm predictions, etc. When the axion-photon coupling is taken into account, 
these results require serious modifications, since in this case the electromagnetic wave 
propagate in chiral (quasi)medium. For instance, in classical terrestrial electrodynamics 
the so-called E- waves and if- waves are introduced [53], for which the radial components 
of the magnetic and electric fields, respectively, are equal to zero. From the mathematical 
point of view, such decomposition of the total electromagnetic field essentially facilitates 
the analysis of the problem. In the axion electrodynamics the introduction of pure E- wave 
or pure H- wave is not possible, since the axionic environment rotates the polarization of 
the electromagnetic waves. 

Our goal is to find some analogs of the Debye potentials for the model, which takes 
into account the axion - photon coupling. Specific feature of our approach is that the 
pseudoscalar field is considered to be formed by relic cosmological axions. This means, 
first, that the axion field <fi is specified, and thus the corresponding evolutionary equa- 
tion for the pseudoscalar field is decoupled from the total set of equations of the axion 
electrodynamics. This approximation is well-motivated, since the density of dark matter 
axions is much higher than the number of axions produced by terrestrial electromagnetic 
field itself. The second specific detail, which also is well-motivated in the context of 
relic cosmological axions, is that the pseudoscalar field (p(t) evolves with time much more 
slowly than the typical terrestrial electric and magnetic fields oscillate. More precisely, in 
our model we consider the time derivative of the pseudoscalar field as a constant, and 
introduce appropriate quantity q, which is directly connected with actual mass density of 
the axionic dark matter in the vicinity of the Earth. 

In other words, we focus on a typical problem of mathematical physics, which came 
from the theory of evolution of the terrestrial electric and magnetic fields, the axion - 
photon coupling being the new key element of theory. 



1.3 Organization of the paper 

The paper is organized as follows. In Section 2 we formulate the model equations of axion 
electrodynamics for the case, when the pseudoscalar field is formed by relic cosmological 
axions born in the early Universe and is revealed now as axionic dark matter. In Subsec- 
tion 2.3 we consider two illustrations: one exact and one approximate solutions to these 
equations, which show that relic axions can form specific electrodynamic configurations 
with parallel coupled electric and magnetic fields. In Section 3 we introduce the so-called 
truncated ((^-independent) model and represent the reduced model equations in spherical 
coordinates. Section 4 is devoted to the analysis of static solutions to the basic equations: 
we discuss the axion electrostatics in Subsection 4.1, the axion magnetostatics in Subsec- 
tion 4.2. In Subsection 4.3 we reduce the obtained results to the dipole-type model of the 
terrestrial magnetic field, and consider some properties of this model in Subsection 4.4. 
In Section 5 we analyze time-dependent solutions to the basic equations of the truncated 
(99-independent) model in the spherical resonator: in Subsection 5.1 we obtain eigen- 
functions for the boundary-value problem, frequencies and amplitudes of the axionically 
coupled modes; in Subsection 5.2 we solve the problem for the case, when the azimuthal 
electric current induces forced oscillations. In Section 6 we focus on the description of 
axionically coupled electric and magnetic oscillations near the Earth surface. Section 7 is 
devoted to discussion of the obtained results. 



2 Axion-photon coupling in the axionic dark matter 
environment 

2.1 The model 

Let us start with the action functional 

S J- J d ^ x ^^^^ + ^ FmnF rnn + ^ KnF rnn_^2 [ Vm 0V"^-V(^)] j (l) 

Here, g is the determinant of the metric tensor g ik , V m is the covariant derivative, R is 
the Ricci scalar, k = G is the gravitational Newtonian coupling constant, A is the 
cosmological constant, % is the Planck constant and c is the speed of light in vacuum. 
The Maxwell tensor F mn is given by 

F mn = V m A n -V n A m , V k F* lk = 0, (2) 

where A m is a potential four-vector of the macroscopic electromagnetic field; F* mn = 
\e mn P q F pq is the tensor dual to F pq ; e mnpq = -^£™ M is the Levi-Civita tensor, E mnpq is 



the absolutely antisymmetric Levi-Civita symbol with _E 0123 =1. 

The first term in (pQ) is the Hilbert-Einstein Lagrangian for the gravity field with 
the account of the cosmo logical constant. The second term \F mn F mn relates to the 
standard linear electrodynamics in vacuum. The term ^<f)F£ in F mn in ([T]) describes the 
pseudoscalar-photon interaction [26]. The symbol stands for an effective pseudoscalar 
field; this quantity is dimensionless providing the second and third terms in the integral 
to be of the same dimensionality. The macroscopic effective axion field, $, is considered 
to be proportional to this quantity $=ty (f>, and the constant ^ is the subject of a special 
discussion. The parameter \I/q has the dimensionality of a force, and its estimation depends 
substantially on the model of the axion coupling, on the value of the axion-photon coupling 
constant gt a 77 ), on the individual axion mass m(A), on the number of axions in the unit 
volume jV(a)) etc. (see, e.g., [TUl EE])- The function V(<p 2 ) describes the potential of the 
pseudoscalar field; in the simplest case we use the quadratic potential V^=/i 2 2 , where 
^_ m ( A)C ^ e dimensionality of inverse length. 

Below we assume that the gravitational background is given, the space-time being of 
the Friedmann - Lemaitre - Robertson - Walker (FLRW) type 

ds 2 = aV) [(dx ) 2 - dl 2 ] , (3) 

with the scale factor a and the Hubble function 

jj ^ da a 

a? dx° a 

Here and below the dot denotes the derivative with respect to time t, which is connected 
with x° by the differential relation 

a(x°) dx° = cdt . (5) 

The metric in the three-space is assumed to be represented in the spherical coordinates 

dl 2 = dr 2 + r 2 {d8 2 + sin 2 9d<p 2 ) , (6) 

and we use the formula y/^g = a 4 r 2 sin6 l . This means that we do not consider back- 
reaction of the electromagnetic field on the gravitational field; moreover, we neglect the 
direct influence of the terrestrial gravity field on the electric and geomagnetic fields in 
comparison with the influence of the relic cosmological axions. 

The effective macroscopic pseudoscalar field, appearing in the action functional ([TJ, 
satisfies the equation 

V fc V fc + = -^2KnF mn , (7) 



however, we assume that the number of axions produced by the macroscopic electromag- 
netic field Fik is much less than the number of relic (primordial) axions created in the 
early Universe. This means that we neglect the electromagnetic source in the right-hand 
side of the equation (j7j) and consider the function <fi(t) to satisfy the decoupled equation 



2„2, 



+ 3H(p + n l c 







The last and very important element of this macroscopic model comprises the relation 
between the effective pseudoscalar field and the dark matter state functions, the energy 
density W(t) and pressure P(t). These relations appear when we put equal the effective 
stress-energy tensor of the macroscopic pseudoscalar field on the one hand, and the phe- 
nomenological stress-energy tensor of the dark matter fluid, on the other hand (see, e.g., 
[37. 251 for details): they have the following form 



Q i 2 jl2 



1 2/2 

a 1 



W . 



P. 



The quantity q is defined as 



and can be readily found as 



dx° 



a 

c 



q 



±—Vw + p. 

For the cold dark matter P — > and W — > pc 2 , thus 

, ca 

q ± -frv/ ? ( DM ) ' 

where P(dm) is the mass density of the cold dark matter. 



(9) 
(10) 

(11) 
(12) 

(13) 



2.2 Electrodynamic equations 

The equations of the axion electrodynamics have the following form 



F ik +(f)F' 



ik 



C 



where J 1 is the electric current four-vector. Taking into account the equations 
can rewrite these equations as follows: 



V k F* = 



J l + J, 



(A) 



where 



—J[ A) = F^ k V^ 



(14) 
one 

(15) 

(16) 



When the electric current four-vector J 1 describes the standard conductivity, we can write 

ji = aE i = a p lk U k = \{g m g kn - g m g km ) a U k F mn , (17) 

where U k is the four- vector of the macroscopic velocity of the system as a whole. Similarly, 
we can write 

— J'iA) = h lkmn V k 0F mn , (18) 

and link this term with axionic conductivity. Clearly, the four-pseudovector -^Vk<P re- 
places the four- vector crU^, when the axionic conductivity is discussed. 



2.3 Prologue: 3-dimensional representation of the model and 
spatially homogeneous solutions 

In 1987 Wilczek (see the paper [27]) considered the application of axion electrodynam- 
ics to the problem of spherically symmetric static dyons. The dyon-type solution was 
obtained using two ingredients: a magnetic monopole plus axion field. In this case the 
static axion field produces static radial electric field using the radial magnetic field of 
the monopole. In fact, Wilczek presented the first (static) example of the longitudinal 
electro-magnetic cluster, since the static electric and magnetic fields are parallel in this 
spherically symmetric configuration. Now we illustrate the idea of longitudinal clusters by 
the example, in which the electric, magnetic and axion fields are spatially homogeneous, 
but depend on time. 

Equations of the axion electrodynamics can be written in the three-vectorial form as 
follows: 

1 d - - - - 

rotE = — —B, divE = 4vrp - B ■ V0 , (19) 
c at 

div£? = 0, rotB = -^-E + — J + -B^-<p-Ex V</>. (20) 
cat c c at 

We assume here that the pseudoscalar (axion) field depends on time only, i.e., <fi=<f)(t), and 
V0=O. Also, we assume that p=0 and J=0, i.e., we study electrodynamic effects in the 
regions out of the electric charges. Then there exists the following exact constant solution: 
the magnetic field is constant, i.e., B=B =const, and the electric field is homogeneous, 
i.e., E=E(t). Indeed, all the equations (|T9|) and (|20|) are satisfied, if 

° = it s + s ^- < 21 > 

Clearly, the exact solution for the electric field reads 

E(t)=E(0)-[<f>(t)-<f>(0))B . (22) 



This means that in the presence of constant magnetic field in addition to the initial arbi- 
trary directed electric field E(0) the electric field collinear to this magnetic field appears 
due to the influence of the evolving system of relic axions. 

Now, let us consider an illustrative approximate model, in which the magnetic field is 
not constant, B=B(t), but it varies with time very slowly. Such situation is typical for the 
problem of infra-low-frequency variations of the terrestrial electric and magnetic fields. 
For instance, the variations of the geomagnetic field on the frequency u/27r=lO~ 3 Hz have 
the effective wavelength of the order A=300 • 10 6 km, which is the distance to the Sun and 
back. Thus, the right-hand side of the first equation in f fl9|) can be considered negligible, 
and the electric field, again, can be treated as homogeneous. Then, one obtains from f l20|) 

E(t) = E(0) -cC drqMBM . (23) 
Jo 

In comparison with the infra-low-frequency variations of the geomagnetic field the function 
q(t) behaves as a constant, and we can use the simplified formula 

B(t) = B + 5B- cos tut -> E(t) = 3(0) - - <p(0)}B - —5B ■ smut, (24) 

CO 

where 6B denotes the amplitude of the geomagnetic field variation, the quantity u> being 
its frequency. This is a good illustration of our idea to focus on the analysis of the 
electric field variations produced by the geomagnetic field variations in the relic axion 
background. Clearly, this electric field induced by the axions is parallel to the geomagnetic 
field variation, that is why we use the term "longitudinal electro-magnetic cluster" to 
indicate this specific mode in the combined electric and magnetic fields variations. 

3 Electrodynamics in the relic axion background 

3.1 Truncated ((^-independent) model 

Let us now apply the model to the case of terrestrial electric and magnetic fields. In fact, 
the geomagnetic field is non-homogeneous, and we have to take into account the curva- 
ture of the magnetic field lines when consider global effects in the Earth's magnetosphere. 
Below we discuss the non-stationary electric and magnetic fields, and start with the poten- 
tials, which do not depend on the azimuthal variable (p, i.e., Ai=Ai(t, r, 9) for i=0, r, 6, if. 
This model is truncated, nevertheless, it illustrates all the important properties of the 
axion-photon coupling in the Earth's electrodynamic system. 

One can introduce the truncated ((^-independent) model using the following approach. 
Let us remind that the flat space-time metric admits the existence of the Killing vector 



£/ \ = 5^ and the Lie derivative of the metric is equal to zero, 

Q Q 

so in the standard spherical coordinates the metric does not depend on (p. Let us assume 
that the potential four- vector inherits this symmetry, i.e., 

= 4)^ + ^4) = °- ( 26 ) 

This provides the condition -^-=0. 
3.2 Master equations 

The equations (114j) can be written as follows. The first equation (for i—0) 

1 d 1 d 

- — (r 2 F 0r ) + -^— — ( S m6F oe ) = 4np, (27) 
r 2 Or r 2 sm 9 09 

4 „ 

where p = — J , does not contain the quantity q. Three equations for i—r 7 9,(p give, 
respectively, 

J^F 0r + - T i— |-(sin £F re ) = -y^— - -r , (28) 
ox u r z sm 6/ atf r z sin ft c 



4- Feo + ^F re = -3-Frv + -r 2 I e , (29) 
ox or smv c 



+ + sin J- ( = g sin 0F r . e - -r 2 /^ sin 2 9 . (30) 

&c u or 09 \r 2 sm9J c 



Here the following notations have been used 

I r = a 4 J r , I e = a 4 J e , I* = a 4 J ip , (31) 

in these terms the equations (j2"T|) -(l3"Ul) do not contain the scale factor a(x°) at all, so that 
the only quantity q reminds that these are the model equations of axion electrodynamics. 

The components of the electric current four-vector are linked by the charge conserva- 
tion law 

c 7ToP + 4|"( r2/r ) + ^|j( sin ^ /9 ) = • ( 32 ) 
Ox" r 2 or sm 9 09 

The second subset of equations in fl2]) gives one equation of the form 

£ F * + ! F » + 5? F " = - (33) 

which does not include the symbol <p. Other equations 

J^V + |:^ = o, j|^ + ^> = o, |:^o + ^^ = o, (34) 

contain only two terms for the truncated model under consideration. 

Linearity of the electrodynamic equations allows us to consider separately static and 
time-dependent solutions of the model. Let us start with the static ones. 



4 Static electric and magnetic fields 



Since the axion field changes extremely slow, and we suppose q=const, it is convenient 
to reconstruct the electric and magnetic fields in terms of two potential components, 
namely, A (r, 6) and A v (r, 6). Let us demonstrate the procedure and the obtained results 
for electric and magnetic fields out of domain of charges distribution. 

4.1 Axion electrostatics 

As usual, we distinguish the coordinate and physical (with subscripts in parentheses) 
components of the electric field, defined as follows 



E T = prO = _0 = p^ ^ ^ 



E e ee = -1^° ee l£ (merid) , (36) 
r z do r 



E 9 = F*> = —^t-^ = —.E^ = . (37) 
r z sin a oip r sin 6 

Let us remind the motivation of these definitions. Since there exist two sets of quantities: 
the covariant and contravariant components of the electric field three- vector (E r , Eg, E v 
and E r ,E 9 ,E lfi , respectively), one needs to introduce the so-called physical components, 
which are assumed to be measured in experiments. We follow the Synge's concept (see, 
e.g., [38J) and introduce the tetrad four- vectors 

A (0) = ~ S > A (rad) = -~K > A (merid) = ' A (azim) = ~~~^Q^ ' ( 38 ) 

which are orthogonal one to another and satisfy the normalization conditions 

fi , ifcA( )A( C ) = l , P«;A(rad)^(rad) = — 1 > ^^(merid) ^(merid) = ~ 1 J 9ik -^(azirn) ^(azim) = — 1 J 

(39) 

when the space-time metric has the form ([3]) with ([6]). The corresponding components of 
the electric field three-vector are defined as 

-E'(rad) = --^j^(rad) > -^(mcrid) = -^i^(merid) ' ^(azim) = ~-^i^(azim) ' (40) 

Cb Cb Hi 

providing the relationships ()35l)-(l371). In terms of standard three-dimensional metric ([6]) 
these quantities can be formally expressed as 



-E(rad) = y\E r E r \ , -^(merid) = \J \E d Eg\ , -E'(azim) = ylE^E^ . (41) 

The equation (1271) with p=0 yields the truncated Laplace equation 

A(o)A = 0. (42) 



Here and below we use the following truncated Laplace operator for the sake of conve- 
nience: 



A 



= f r 2 — ) + — 



1 <9 / . o> \ m 2 
sln~#d# V d# j ~ sin 2 # 



(43) 



The equation (1421) does not contain the quantity q, so, the corresponding solution is 
well-known in classical mathematical physics. E.g., for the domain r > R it has the form 



n=0 ^ r ' 



(44) 



where R is the Earth radius. We use the standard definitions: P n (x) for the Legendre 
polynomials [39j EH] and P^ m ^(x) = (1— x 2 ) 2 ( -P ra (x) for the adjoint Legendre poly- 
nomials. As usual, the coefficients C n can be found from boundary conditions (we do not 
consider this problem here). Clearly, the static electric field components 



E 



(rad) 



E 

n=0 



[11 



^C n (f) P„(cos ( 



E 



(merid) 



Si ^ 



n=0 



COS I 



(azim) 



do not feel the influence of axions. 



(45) 
(46) 



4.2 Axion magnet ost at ics 

Magnetic field components are defined as follows: 

1 



1 9A V _ B 
r 2 sin# 90 , " 1 



B° = E 



r 2 sin 6 
1 



r rl „ , -D (merid) • /i o ' 

r sm 6 or 



dA e dA r 



r 2 sin 6 



E H 



B 



(azim) 



dr 86 



where the so-called physical components of the magnetic field 



P(rad) = ^\B r B r 



B 



(merid) 



BaB b 



(azim) 



B v B<p\ 



are obtained in analogy with the electric field components (see (155 ]) - (j4"l~jT 
equations (|28|) and ( )29|) convert now into 

— (sin 6F rd - qA v ) = , — (sin 6F rd - qA v ) = , 



we obtain immediately that 



g 4 

sine/ 



(47) 

(48) 
(49) 

(50) 
Since the 

(51) 
(52) 



Then the equation (130)) reduces to 











and using the new potential U(r, 9), defined by 



A w = r sin 9 U , U 



{ 1 8 ) 


2 2 


\sin9d9 J 


+ q r 



A ip = 0, 



(53) 



A v Af\ 



we obtain the truncated ((^-independent) Helmholtz equation 

A {1) U + q 2 U = 0. 

The solution of (|55p has the form 

oo 

f/=£P«(cos#) 5R n (r,g), 



(54) 



(55) 



(56) 



n=l 



where the radial function 9ft n (r, g) depends essentially on the parameter q. When q = 0, 
i.e., we neglect the axion influence at all, this function is [39| 135] 

r \ — (ra+l) 

(57) 



When q ^ 0, the radial function can be expressed in terms of Bessel functions 



^n(r,q) = \j— {An 



J n+ dqr) + 



71 



-qR 



J -(n + k)(l r ) 



r(n+|) V2 

where J n+ i{qr) is the Bessel function of the first kind with the half-integer index 
The constants of integration are chosen so that lim^o 9ft„(r, q) = 9ft n (r, 0) (see (1ST)) ). 
The corresponding physical components of the magnetic field are 



(5f 



5, 



(rad) 



n=l 



1 d 



B (melid) = -^P«(cos0)--[rl? n (r, g )] 

n=l r ^ r 



5 (azim ) = -g E P i 1} ( cos ^)^n(r, q) 

n=l 



(59) 

(60) 
(61) 



The main new feature is the following: the azimuthal component -B( az im), being equal to 
zero at q=0, becomes nonvanishing at q ^ 0. 



4.3 Dipole-type magnetic field 

In the zeroth-order approximation the terrestrial magnetic field can be described in terms 
of geomagnetic dipole. In the framework of this model we obtain basic formulas if put 
n = 1 into (!56j) - fl6Tl) . what yields 

2 

5( rad ) = -- cos ® ^i( r ,Q) , (62) 



^(merid) = ~~ [ r ^l( r ,?)] , (63) 



sin 6* d 

r dr 

£( az im) = -q sin 9 3?i (r, q) . (64) 



We consider two auxiliary angles. The first auxiliary angle 7, defined as 

tan 7 = ^ azim) /p? = ^tang, (65) 

is linear in the parameter q and depends on the meridional angle 9. As for the second 
auxiliary angle 5, defined as 

, 1 -1 



tan S = L T ^ = QR 



1 | p gjWg) 



(66) 



and thus 



it does not depend on 9. For the domain r > Rby substitutions A n = and B\ = we 
obtain 

( r ) — ^ ( cos <? r + <? r s * n 9 r ) > (67) 

-B(rad) ( r ? 9) = — j cos 9 (cos gr + qr sin gr) , (68) 

/isin^ r 2 2 ] 

-D(merid) {r, q) = 3 — (cos gr + qr sin gr) — g r cos qr , (69) 

5( az im) ( r ' q) = -q sin ( cos <? r + <? r sin <? r ) • (70) 
In the approximation g — >■ these formulas give 

■B(rad)(r, 0) = -~ cos6» , 5 (merid )(r, 0) = /i ^ 6> , B {azim) (r, 0) = , (71) 

as it should be in classical theory of static dipole geomagnetic field. Let us emphasize 
that in the approximation linear in q the angle 5 is 

tan 5 = -qR, (72) 



thus, the axion-photon coupling can be considered as one of the explanations of the 
(visual) magnetic poles drift. 



4.4 Axionically coupled magnetic field as a function of altitude 

Classical dipole-type magnetic field decreases monotonically at r — » oo according to the 
law 1/r 3 (see, e.g., (I7T1)). In the axionic environment, i.e., when q ^ 0, the magnetic 
field components behave non- monotonically (see (ITU]) ). Indeed, there is infinite number of 
values of the radius r=RJ m -., for which the radial function 3?i(r, q) takes zero value; they 
can be found from the equation 

cos [qR* m) ] + qR\ m) sin [qR\ m) ] = , (73) 

or equivalent ly, 

ctgz = -z , z = qR\ m) . (74) 

At the same altitudes r=R*r m \ the radial and azimuthal components of the magnetic field 
change the signs. The position of the first null is defined by the inequalities z < ti 
and R*^ < ir/q. The radial function has infinite number of extrema at r = RTfy, these 
quantities satisfy the equation 

tgz = ^-- z , z = qR* { ; ) . (75) 

The first extremum is placed at z < 3tt/2, Rt^ < 3n/2q. In other words, the relic axions 
change essentially the configuration of the magnetic field. 

5 Time dependent electric and magnetic fields in the 
relic axion background 

Taking into account the linearity of the model of terrestrial axion electrodynamics in the 
relic axion background, we distinguish two cases in the non-stationary problem. First 
of all, we analyze natural oscillations and consider the electric current four-vector to be 
vanishing. Then we discuss the effects induced by the azimuthal electric current. 

5.1 Natural oscillations of electric and magnetic fields in the 
spherical resonator Earth-Ionosphere 

5.1.1 Key equations 

When the electric current four- vector is equal to zero, the equations (I27l)-(l30l) and ( 1331) . 
(jMj) can be reduced to the coupled pair of equations for two new potentials U and V, 



which play in axion electrodynamics the same role as the Debye potentials [35] play in 
the electrodynamics of Faraday - Maxwell. Let us use the substitutions 

For = ^ jWsintf) = £ (rad) , (76) 

r sin d ad 



Feo = ~fr( rV ) = - rE (merid) , (77) 

d 

F v0 = ~r sin 9— U = -r sm 6>£ (azim) , (78) 



Fro = r[qU + ^-Vj = -r£ (azim) , (79) 
d 

F r!fi = sin 6—(rU) = -r sin 9B {merid) , (80) 

F Q9 = r-^(U sin6) = -r 2 sm6B {r&d) , (81) 
A (p = rUshi6. (82) 
Then the equations (127)1 . (1251) . (129"]) and (13"lj) turn into identities, the equation (13"3"]) gives 

^»"-^ = ^. 
and the equations (130]) takes the form 

^"^ = -^-^- (84) 

We indicate the pair of equations (153]) and (1841) as key equations. When q ^ 0, using the 
decoupling procedure for the equations (183]) and ( 184]) . we obtain the key equation of the 
fourth order for the function V: 



A (D- 



d 2 12 



+ g 2 A (1) ^ = 0. (85) 



dx 02 

When V is found, the function U can be extracted from ( 183]) by integration over time. 



5.1.2 Boundary value problem, eigenvalues and eigenfunctions 

Let us consider electromagnetic oscillations in the spherical resonator bounded by the 
Earth surface (r=R) and the bottom edge of the Earth Ionosphere (r=R if ) in the frame- 
work of the truncated model. 

Let the quantities jF n j(r,9) be the eigenfunctions of the operator A^) satisfying the 
boundary conditions 

F nj {R,6) = 0, J c nj {R*,6) = 0. (86) 



They have the well-known multiplicative form 



J^(r,0) = P«(cos0) H nj (r), 



»7) 



where the radial functions 7i n j(r) are 



J n + I J -(n + |) (^) - J n + | (^) -L^) 



(88) 



Boundary conditions f lSBj) are satisfied, when the parameters are found from the 
equation 



Jn+l {vfR) JL M ) (^*) = (^*) ^(n + |) ("j"'*) 



39) 



where the index j = 0, 1, 2, ... counts the positive zeros of the equations ([BSD - Of course, 
the boundary conditions (|86|) are chosen as the simplest example, which can illustrate 
our main statements; as for applications, we intend to use more sophisticated boundary 
conditions in special papers. 

We search for the potential functions U and V in the form 



U(t,r,e) = EES'W ^ (<*»*) H nj (r) 

n=0j=0 



(90) 



V(i,r,e) = EE^W ^(cosfl) W fti (r) 

n=0 j=0 



(91) 



where we put for the sake of simplicity that x° = ct. The key equations (l84"j) and ( 183)) 
give the following coupled pair of equations for the time dependent amplitudes: 



f OA 2 2 



u nj - qcv nj = 



v nj + c 2 (vj n) ) v nj + qcu nj = . 
The axionic parameter enters these equation non-equivalently. 



(92) 
(93) 



5.1.3 Solutions for the mode amplitudes 

In order to simplify the analysis we suppose here that the axionic parameter q is positive, 
q > 0. The system of coupled equations (1921) and (1931) has standardly the solutions of three 
types. When q is so small that q < za- for all j and n, we obtain harmonic oscillations 
for both functions 



v n j{t) 



C[Jj cos uJi n jt + C^Jj sin LUi n jt 



+ 



C lr lj COS UJ2njt + C^nj SID. UJ2njt 



(94) 



and 



u nj (i)-- 



1 9 

1 -tjT) 



CoJi cosu} lnit — Cl': Sill Uinjt 



(1) 



+ 



1 _i 2_ 

it (n) 



Cj n j Sin 0J2 n jt — C2nj COS W2raj^ 
.(1) rrW ^(2) OT1 i ^(i) 



(95) 



The constants C) nj -, C^-, Ci n j and can be found from initial data. The hybrid 
frequencies 

(96) 



(n) 



\ 



V 



(n) 



^2nj — cv j 



(n) 



\ 



1 + 



(n) 



coincide when g=0, i.e., when there is no axion-photon coupling. 

Formally speaking, for some mode numbers, say, n < and j < j*, the quantity q 
can exceed v!j™ . The corresponding modes are non-harmonic and one should replace the 
first pair cosui n ji and smoji n jt by hyperbolic functions coshr ln ,,t and sinhr lry t, where 

r ln ,- = c ^ n) ./^-i. 

,(n.) 



When q—u^* , the corresponding mode amplitudes 



u*(0)- 



u*(0) 



+ 



cq 

«*(*) -- 
1 



+ 



v*(0)- 



u*(0) 



eg 



y — ~ V*(0) I- ~ 

cos v 2cgt H — p — sin v 2cgt . 



2^2 



m*(0)-1 
w*(0)4 



u*(0) 



2cg 

«*(0)~ 



cgt 
IT 



u*(0)- 



\/2cg 



(97) 



eg 



+ 



eg 



sin y/2cqt — * cos \Plcqt , 
2 eg 



(98) 



display the possibility of oscillations on the specific (axionic) frequency WA=V / 2cg, as well 
as, the linear growth with time of the [/-potential. 



5.2 Forced oscillations produced by an azimuthal electric cur- 
rent 

Let us consider the models in which the only azimuthal component of the electric current 
four- vector J v is no n- vanishing. This case is very illustrative just for the truncated ((p- 
independent) model. Indeed, the charge conservation law (132]) is now satisfied identically. 
The electrodynamic equations (1551) . (154]) with the definitions (176"]) - (152]) transform into 



d 2 d 
d 2 d 4tt 



(99) 
(100) 



where 



/ = J\IfL 



I^r sin 6 . 



Decoupled equation for V is modified accordingly as 



2 



dx 02 



i) 



c dx° 



(101) 



(102) 



Repeating the decomposition of the potentials V and U fulfilled in Subsection 5.1., we 
obtain the differential equations 

, 2 



Q 



u 



qcv nj = I nj (t) 



v nj + c 2 (vj n) ) v nj + qcii nj = , 
where the functions I n j(t) denote the coefficients in the decomposition 



47T r .~ 

—I(t,r, 
c 



(103) 
(104) 

(105) 



= EE J *(*) p £ ) («»0)^ 

n=0 j=0 

When q— 0, the key equations are not coupled, so the only U potential feels the influence 
of the stimulating current of this type. The case q ^ is much more sophisticated. 

The particular solutions v n j(t), which satisfy fl 1 3 [) and (1104p and null initial data 
v nj{®) = = v n j(0), have the following form 



V2Juj 2 +u 2 



dTl nj (r) cosu 2n j(t-T) - cos u lnj (t-T) 

J lnj^ rUJ 2nj ""' 

Clearly, this solution vanishes at q — > 0, since in this limit U2 n j ^i n j- 

The particular solutions u* nj (t) with null initial data u n A0) = = u n AQ) is 



(106) 



1 ft 



2 Jo 



smcui n j(t—T) smui2nj(t—T~) 



(107) 

Thus, when q=0, the azimuthal current produces stimulated azimuthal electric field, 
meridional and radial magnetic field. When q ^ 0, additional radial and meridional 
electric field, as well as the azimuthal magnetic field also appears as a result of the axion- 
photon coupling. When the azimuthal component of the electric current is periodic, and 
the frequency is, say, Q, i.e., 

hjii) = Inj sin fit, (108) 

we obtain stimulated oscillations on the heterodyne frequencies U\ n j ± Q and Wmj ± fi- 
For instance, for the amplitude u^At) one obtains: 



4:UJ 2 nj U% nj 



I n j 1 2 (wi nj + w 2nj) sin ^t+^lnj sin (^i nj -fi)t- sin (ui nj +Q)t 



+00 



Inj 



sin (uj2nj— tyt — sin (oj2 n j+^)t } 



+ 
(109) 



For the mode amplitude v n j(t) the result is similar. 



6 An application: axionically coupled electric and 
magnetic fields near the Earth's surface 



6.1 Longitudinal cluster formed by meridional electric and mag- 
netic fields 

Let us consider the electric and magnetic fields at the Earth surface r = R, using the 
definitions ( |76|) - (l8T!) and the solutions (!90|) - (l96|) to the key equations. The radial and 
azimuthal components of the electric field in this model vanish at the Earth surface 

-E'(rad) (£, -R, 0) = , -E(azim) (t, -R, 9) = . (HO) 

The radial and azimuthal components of the magnetic field also take zero values 

B (nd )(i, R,9) = 0, B [azim) {t, R,9) = 0. (Ill) 
Only meridional components of electric field 

oo oo 

E (merid) (t,R,9) = P«(cos0) g nj (R) , (112) 

71=0 j=0 

and magnetic field 

oo oo 

B {mcrid) (t,R,9) = -J2Y, u ni(t) Pi 1) (™s9) g nj {R), (113) 

n=0 j=0 

are non-vanishing and thus attract attention in this model. In the formulas (11121) and 
(11131) we used the auxiliary function 



GnAR)=% J' nH & n) R) J. M) ti n) R) -J nH ^R) f_ M) {vfR 



(114) 



Clearly, the meridional electric and magnetic fields are coupled by the axion field, and 
we deal with longitudinal cluster of electromagnetic modes. Let us analyze the time 
dependent variations of these fields by two examples. 

6.2 Perturbations of meridional electric field produce variations 
of the meridional magnetic field 

Let us consider, first, the situation when initial data are the following 

v ni (0)^0, t> ni (0)^0, 1^(0) =0, u nj {0) = 0. (115) 



This means that initially only the meridional electric field is perturbed. The corresponding 
solution for the v n j amplitude is (see fl9] 



V n j(t) = 2 U nj(0) COS Ui nj t + COS U 2n jt 



+ 



+ ^nj(O) — 77 7(n)j [wiry- sin u;i ni t + c; 2ni smu 2n jt] 



(116) 



2 " ' ' UJlnj<^2nj 

The solution for the M nj - amplitude gives now no n- vanishing result, though these functions 
started from zero initial values (see (I115P ): 

1 



u nj (t) 



2y/Ul n jU2nj 



7(„)j {vnj (0) ui n j sin uj 2n jt - io 2n j sin U lnj t 



+ 



+i>nj(0hf n )j COS Uinjt - COS UJ 2n jt } • 



Here we introduced the following convenient constants: 



7(n)j 



on 



-2 



;ii7) 



:ii8i 



For the resonant case, when q=v 1 j* for some mode numbers n* and j*, one can illustrate 
the t/-mode amplitude generation by the formula 



u*(t) 



2V2 



sin \plcqt — v2cqi 



u„(0) . 2 cqt 

H sm — = . 

cq ^2 



(119) 



Clearly, there is the term linear in time, and for small value of cqt the expression fll 19j) 
has the leading order term oc cqt 2 . 

If q = 0, so that c<;i n j=c<;2nj = ^o n j, one obtains that u n j(t) = 0, thus the magnetic 
field remains unperturbed. In the approximation linear with respect to q one obtains 



u nj (t) 



qc 



2co, 



Onj 



Vnj(t) = fnj(0) COS UlQ n jt + SIB. U 0nj t , 

[v n j(0) OJQnjt COS UJ 0nj t - sin UjQnjt + V n j(0)t SU1 Uonjt} ■ 



(120) 



(121) 



In particular, we obtain very convenient formula for estimations 



U nj{t) = -QCt V nj (t) 



:i22) 



when the initial data have specific form v n j(0) = 0. 



6.3 Perturbations of meridional magnetic field produce varia- 
tions of the meridional electric field 

Let us consider now the situation when initial data are the following 



u nj (0)^0, u nj (0)^0, v nj (0) = 0, v nj (0) = 0. 



(123) 



This means that initially only the meridional magnetic field is perturbed. Now the v n j 
amplitudes become non-vanishing at t > due to the axion-photon coupling: 



v nj (t) 



7(n)j 



, \u nj (0) COS U 2nj t - COS LU lnj t 

^^&l n jU! 2n j L L 

+Unj(0)lf n )j [u}2njSmU lnj t - Ui nj 8mu 2n jt } • 



+ 



(124) 



The u n j amplitudes evolve as follows: 



u nj (t) 



1 



{u n j{0)lf n)j [u>l nj COS U lnj t + U 2 lnj COS U 2nj t 



2Ui n jU2nj 

+u nj (0) \u 2nj sin Uinjt + u lnj sin uj 2n jt] } . 



+ 



(125) 



Again, when g=z/™*, one can illustrate the f-mode amplitude generation by the formulas 



,r\ ^*(0) . 2 cqt «*(0) 



sin \/2cqt + V2cqt 



(126) 



The term linear in time appears now in u*(t). 

Clearly, if q=0, so that Ui n j=u 2n j, one obtains that v n j(t)—0, thus the electric field 
remains unperturbed. Analogously, in the linear approximation we obtain 



u nj (t) = u nj (0) cosuj 0nj t + Vn3 ° sin u 0nj t , 



(127) 



v nj (t) 



qc 



2u>, 



Onj 



with 



{u nj (0) UjQnjt COS Uonjt ~ sin Wonji 



u nj(*) = ~2 qCt Un ^ ' 



+ M nj (0)tsina;onii} , (128) 

(129) 



when the initial data are chosen as u n j(0) = 0. 



7 Discussion 

1. We studied one model of coupling of terrestrial magnetic and electric fields with a relic 
axion background. Exact solutions obtained in the framework of the axion electrodynam- 
ics illustrate three important conclusions. 



1.1. Relic dark matter axions produce in the terrestrial electrodynamic system oscilla- 
tions of a new type, which belong to the class of Longitudinal Electro-Magnetic Clusters 
(LEMCs). LEMC has the following specific feature: the electric and magnetic fields are 
parallel to one another and are coupled by axionic field (see, e.g., the formulas (jllOp - 
(11141) ); in the absence of axions such oscillations decouple. Electric and magnetic fields 
of this type are correlated: the oscillations are characterized by identical frequencies but 
different phases (see, e.g., the formulas (II 16ft - (II 171) ) and (I124p - (ll25p )). Generally, there 
are two sets of hybrid frequencies of LEMCs (see the formulas (M)|) ). and these frequencies 
coincide when the axionic dark matter influence is negligible. 

1.2. Relic dark matter axions deform the static terrestrial magnetic field. First, 
in case when the original geomagnetic field has the radial and meridional components 
only, the axion-photon coupling produces a supplementary azimuthal component (see, 
e.g., (l6Tj) . (170|) ); this effect contributes to the phenomenon of the Earth's magnetic pole 
drift. Second, due to the axion-photon coupling the dependence of the magnetic field 
on the altitude becomes non-monotonic: the surfaces appear, on which the radial and/or 
meridional components of the geomagnetic field change the signs, as well as reach maxima 
and minima (see Subsection I4.4p . 

1.3. The truncated (^-independent) model can be effectively described in terms of two 
electromagnetic potentials U and V, which can be considered as axionic generalizations 
of the well-known Debye potentials. 

2. Estimations of the effect produced by relic axion background in the terrestrial 
electrodynamic system depend on how small the dimensionless parameter £ = qR is. It 
can be estimated as follows. Let us remind that in the laboratory frame of reference 
\q\=i§^^P(dm), and use the natural units, in which h—1 and c— 1. According to [ID] 
the mass density of the dark matter in the Solar system is estimated to be P(dm) — 
0.033 M( Sun )pc~ 3 or in the natural units P(dm) — 1-25 GeV ■ cm -3 . According to [TOl [36] 
the parameter \J/ is reciprocal to the axion-photon-photon coupling constant pA-y-y, i-e., 
^— =pA77 and PA77 itself can be (optimistically) estimated as PA77 — 

lO^GeV -1 . Thus, we 

obtain £ ~ 10~ 7 . This means that the effective frequency of LEMCs is about va = c ( 1 = % — 
0.5 • 10 _5 Hz. This frequency belongs to the range of infra-low frequencies; variations of 
electric and magnetic fields of this type are studying in the experimental group of Vladimir 
University during more than 40 years [HI W2\ |4"3"] . 

3. We consider this paper as theoretical grounds for the program of testing of the 
LEMCs in the Earth's atmosphere. Detailed technical description of the experiments and 
first experimental results are planned to be discussed in a special paper. 
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